A Human Homozygous HELQ Missense Variant Does Not Cause Premature Ovarian Insufficiency in a Mouse Model.

Disruption of meiosis and DNA repair genes is associated with female fertility disorders like premature ovarian insufficiency (POI). In this study, we identified a homozygous missense variant in the HELQ gene (c.596 A>C; p.Gln199Pro) through whole exome sequencing in a POI patient, a condition associated with disrupted ovarian function and female infertility. HELQ, an enzyme involved in DNA repair, plays a crucial role in repairing DNA cross-links and has been linked to germ cell maintenance, fertility, and tumour suppression in mice. To explore the potential association of the HELQ variant with POI, we used CRISPR/Cas9 to create a knock-in mouse model harbouring the equivalent of the human HELQ variant identified in the POI patient. Surprisingly, Helq knock-in mice showed no discernible phenotype, with fertility levels, histological features, and follicle development similar to wild-type mice. Despite the lack of observable effects in mice, the potential role of HELQ in human fertility, especially in the context of POI, should not be dismissed. Larger studies encompassing diverse ethnic populations and alternative functional approaches will be necessary to further examine the role of HELQ in POI. Our results underscore the potential uncertainties associated with genomic variants and the limitations of in vivo animal modelling.

Deficiency of the mitochondrial ribosomal subunit, MRPL50, causes autosomal recessive syndromic premature ovarian insufficiency.

Premature ovarian insufficiency (POI) is a common cause of infertility in women, characterised by amenorrhea and elevated FSH under the age of 40 years. In some cases, POI is syndromic in association with other features such as sensorineural hearing loss in Perrault syndrome. POI is a heterogeneous disease with over 80 causative genes known so far; however, these explain only a minority of cases. Using whole-exome sequencing (WES), we identified a MRPL50 homozygous missense variant (c.335T > A; p.Val112Asp) shared by twin sisters presenting with POI, bilateral high-frequency sensorineural hearing loss, kidney and heart dysfunction. MRPL50 encodes a component of the large subunit of the mitochondrial ribosome. Using quantitative proteomics and western blot analysis on patient fibroblasts, we demonstrated a loss of MRPL50 protein and an associated destabilisation of the large subunit of the mitochondrial ribosome whilst the small subunit was preserved. The mitochondrial ribosome is responsible for the translation of subunits of the mitochondrial oxidative phosphorylation machinery, and we found patient fibroblasts have a mild but significant decrease in the abundance of mitochondrial complex I. These data support a biochemical phenotype associated with MRPL50 variants. We validated the association of MRPL50 with the clinical phenotype by knockdown/knockout of mRpL50 in Drosophila, which resulted abnormal ovarian development. In conclusion, we have shown that a MRPL50 missense variant destabilises the mitochondrial ribosome, leading to oxidative phosphorylation deficiency and syndromic POI, highlighting the importance of mitochondrial support in ovarian development and function.

Premature Ovarian Insufficiency in CLPB Deficiency: Transcriptomic, Proteomic and Phenotypic Insights.

CONTEXT: Premature ovarian insufficiency (POI) is a common form of female infertility that usually presents as an isolated condition but can be part of various genetic syndromes. Early diagnosis and treatment of POI can minimize comorbidity and improve health outcomes. OBJECTIVE: We aimed to determine the genetic cause of syndromic POI, intellectual disability, neutropenia, and cataracts. METHODS: We performed whole-exome sequencing (WES) followed by functional validation via RT-PCR, RNAseq, and quantitative proteomics, as well as clinical update of previously reported patients with variants in the caseinolytic peptidase B (CLPB) gene. RESULTS: We identified causative variants in CLPB, encoding a mitochondrial disaggregase. Variants in this gene are known to cause an autosomal recessive syndrome involving 3-methylglutaconic aciduria, neurological dysfunction, cataracts, and neutropenia that is often fatal in childhood; however, there is likely a reporting bias toward severe cases. Using RNAseq and quantitative proteomics we validated causation and gained insight into genotype:phenotype correlation. Clinical follow-up of patients with CLPB deficiency who survived to adulthood identified POI and infertility as a common postpubertal ailment. CONCLUSION: A novel splicing variant is associated with CLPB deficiency in an individual who survived to adulthood. POI is a common feature of postpubertal female individuals with CLPB deficiency. Patients with CLPB deficiency should be referred to pediatric gynecologists/endocrinologists for prompt POI diagnosis and hormone replacement therapy to minimize associated comorbidities.

SOX Genes and Their Role in Disorders of Sex Development.

SOX genesare master regulatory genes controlling development and are fundamental to the establishment of sex determination in a multitude of organisms. The discovery of the master sex-determining gene SRY in 1990 was pivotal for the understanding of how testis development is initiated in mammals. With this discovery, an entire family of SOX factors were uncovered that play crucial roles in cell fate decisions during development. The importance of SOX genes in human reproductive development is evident from the various disorders of sex development (DSD) upon loss or overexpression of SOX gene function. Here, we review the roles that SOX genes play in gonad development and their involvement in DSD. We start with an overview of sex determination and differentiation, DSDs, and the SOX gene family and function. We then provide detailed information and discussion on SOX genes that have been implicated in DSDs, both at the gene and regulatory level. These include SRY, SOX9, SOX3, SOX8, and SOX10. This review provides insights on the crucial balance of SOX gene expression levels needed for gonad development and maintenance and how changes in these levels can lead to DSDs.

Two ovarian candidate enhancers, identified by time series enhancer RNA analyses, harbor rare genetic variations identified in ovarian insufficiency.

The genetic regulation of ovarian development remains largely unclear. Indeed, in most cases of impaired ovarian development-such as 46,XX disorders of sex development (DSD) without SRY, and premature ovarian insufficiency (POI)-the genetic causes have not been identified, and the vast majority of disease-associated sequence variants could lie within non-coding regulatory sequences. In this study, we aimed to identify enhancers of five ovarian genes known to play key roles in early ovarian development, basing our analysis on the expression of enhancer derived transcripts (eRNAs), which are considered to characterize active enhancers. Temporal expression profile changes in mouse WT1-positive ovarian cells were obtained from cap analysis of gene expression at E13.5, E16.5 and P0. We compared the chronological expression profiles of ovarian-specific eRNA with expression profiles for each of the ovarian-specific genes, yielding two candidate sequences for enhancers of Wnt4 and Rspo1. Both sequences are conserved between mouse and human, and we confirmed their enhancer activities using transient expression assays in murine granulosa cells. Furthermore, by sequencing the region in patients with impaired ovarian development in 24 patients, such as POI, gonadal dysgenesis and 46,XX DSD, we identified rare single nucleotide variants in both sequences. Our results demonstrate that combined analysis of the temporal expression profiles of eRNA and mRNA of target genes presents a powerful tool for locating cis-element enhancers, and a means of identifying disease-associated sequence variants that lie within non-coding regulatory sequences, thus advancing an important unmet need in forward human genetics.

Whole exome sequencing reveals copy number variants in individuals with disorders of sex development.

Complete androgen insensitivity syndrome (CAIS), where 46,XY individuals present as female, is caused by variants in the androgen receptor gene (AR). We analyzed the DNA of a patient with suspected CAIS using a targeted gene sequencing panel and whole exome sequencing (WES) but did not detect any small nucleotide variants in AR. Analysis of WES data using our bioinformatics pipeline designed to detect copy number variations (CNV) uncovered a rare duplication of exon 2 of AR. Using array comparative genomic hybridization, the duplication was found to span 43.6 kb and is predicted to cause a frameshift and loss of AR protein. We confirmed the power of our WES-CNV detection protocol by identifying pathogenic CNVs in FSHR and NR5A1 in previously undiagnosed patients with disorders of sex development. Our findings illustrate the usefulness of CNV analysis in WES data to detect pathogenic genomic changes that may go undetected using only standard analysis protocols.

Generation and mutational analysis of a transgenic mouse model of human SRY.

SRY is the Y-chromosomal gene that determines male sex development in humans and most other mammals. After three decades of study, we still lack a detailed understanding of which domains of the SRY protein are required to engage the pathway of gene activity leading to testis development. Some insight has been gained from the study of genetic variations underlying differences/disorders of sex determination (DSD), but the lack of a system of experimentally generating SRY mutations and studying their consequences in vivo has limited progress in the field. To address this issue, we generated a mouse model carrying a human SRY transgene able to drive testis determination in XX mice. Using CRISPR-Cas9 gene editing, we generated novel genetic modifications in each of SRY's three domains (N-terminal, HMG box, and C-terminal) and performed a detailed analysis of their molecular and cellular effects on embryonic testis development. Our results provide new functional insights unique to human SRY and present a versatile and powerful system in which to functionally analyze variations of SRY including known and novel pathogenic variants found in DSD.

Meiotic genes in premature ovarian insufficiency: variants in HROB and REC8 as likely genetic causes.

Premature ovarian insufficiency (POI), affecting 1 in 100 women, is characterised by loss of ovarian function associated with elevated gonadotropin, before the age of 40. In addition to infertility, patients face increased risk of comorbidities such as heart disease, osteoporosis, cancer and/or early mortality. We used whole exome sequencing to identify the genetic cause of POI in seven women. Each had biallelic candidate variants in genes with a primary role in DNA damage repair and/or meiosis. This includes two genes, REC8 and HROB, not previously associated with autosomal recessive POI. REC8 encodes a component of the cohesin complex and HROB encodes a factor that recruits MCM8/9 for DNA damage repair. In silico analyses, combined with concordant mouse model phenotypes support these as new genetic causes of POI. We also identified novel variants in MCM8, NUP107, STAG3 and HFM1 and a known variant in POF1B. Our study highlights the pivotal role of meiosis in ovarian function. We identify novel variants, consolidate the pathogenicity of variants previously considered of unknown significance, and propose HROB and REC8 variants as new genetic causes while exploring their link to pathogenesis.

Targeting the Non-Coding Genome for the Diagnosis of Disorders of Sex Development.

Disorders of sex development (DSD) are a complex group of conditions with highly variable clinical phenotypes, most often caused by failure of gonadal development. DSD are estimated to occur in around 1.7% of all live births. Whilst the understanding of genes involved in gonad development has increased exponentially, approximately 50% of patients with a DSD remain without a genetic diagnosis, possibly implicating non-coding genomic regions instead. Here, we review how variants in the non-coding genome of DSD patients can be identified using techniques such as array comparative genomic hybridization (CGH) to detect copy number variants (CNVs), and more recently, whole genome sequencing (WGS). Once a CNV in a patient's non-coding genome is identified, putative regulatory elements such as enhancers need to be determined within these vast genomic regions. We will review the available online tools and databases that can be used to refine regions with potential enhancer activity based on chromosomal accessibility, histone modifications, transcription factor binding site analysis, chromatin conformation, and disease association. We will also review the current in vitro and in vivo techniques available to demonstrate the functionality of the identified enhancers. The review concludes with a clinical update on the enhancers linked to DSD.

New insights into the genetic basis of premature ovarian insufficiency: Novel causative variants and candidate genes revealed by genomic sequencing.

Ovarian deficiency, including premature ovarian insufficiency (POI) and diminished ovarian reserve (DOR), represents one of the main causes of female infertility. POI is a genetically heterogeneous condition but current understanding of its genetic basis is far from complete, with the cause remaining unknown in the majority of patients. The genes that regulate DOR have been reported but the genetic basis of DOR has not been explored in depth. Both conditions are likely to lie along a continuum of degrees of decrease in ovarian reserve. We performed genomic analysis via whole exome sequencing (WES) followed by in silico analyses and functional experiments to investigate the genetic cause of ovarian deficiency in ten affected women. We achieved diagnoses for three of them, including the identification of novel variants in STAG3, GDF9, and FANCM. We identified potentially causative FSHR variants in another patient. This is the second report of biallelic GDF9 and FANCM variants, and, combined with functional support, validates these genes as bone fide autosomal recessive "POI genes". We also identified new candidate genes, NRIP1, XPO1, and MACF1. These genes have been linked to ovarian function in mouse, pig, and zebrafish respectively, but never in humans. In the case of NRIP1, we provide functional support for the deleterious nature of the variant via SUMOylation and luciferase/ß-galactosidase reporter assays. Our study provides multiple insights into the genetic basis of POI/DOR. We have further elucidated the involvement of GDF9, FANCM, STAG3 and FSHR in POI pathogenesis, and propose new candidate genes, NRIP1, XPO1, and MACF1, which should be the focus of future studies.

STAG3 homozygous missense variant causes primary ovarian insufficiency and male non-obstructive azoospermia.

Infertility, a global problem affecting up to 15% of couples, can have varied causes ranging from natural ageing to the pathological development or function of the reproductive organs. One form of female infertility is premature ovarian insufficiency (POI), affecting up to 1 in 100 women and characterised by amenorrhoea and elevated FSH before the age of 40. POI can have a genetic basis, with over 50 causative genes identified. Non-obstructive azoospermia (NOA), a form of male infertility characterised by the absence of sperm in semen, has an incidence of 1% and is similarly heterogeneous. The genetic basis of male and female infertility is poorly understood with the majority of cases having no known cause. Here, we study a case of familial infertility including a proband with POI and her brother with NOA. We performed whole-exome sequencing (WES) and identified a homozygous STAG3 missense variant that segregated with infertility. STAG3 encodes a component of the meiosis cohesin complex required for sister chromatid separation. We report the first pathogenic homozygous missense variant in STAG3 and the first STAG3 variant associated with both male and female infertility. We also demonstrate limitations of WES for the analysis of homologous DNA sequences, with this variant being ambiguous or missed by independent WES protocols and its homozygosity only being established via long-range nested PCR.

Analysis of NR5A1 in 142 patients with premature ovarian insufficiency, diminished ovarian reserve, or unexplained infertility.

Ovarian deficiency, including diminished ovarian reserve and premature ovarian insufficiency, represents one of the main causes of female infertility. Little is known of the genetic basis of diminished ovarian reserve, while premature ovarian insufficiency often has a genetic basis, with genes affecting various processes. NR5A1 is a key gene required for gonadal function, and variants are associated with a wide phenotypic spectrum of disorders of sexual development, and are found in 0.26-8% of patients with premature ovarian insufficiency. As there is some debate about the extent of involvement of NR5A1 in the pathogenesis of ovarian deficiency, we performed an in-depth analysis of NR5A1 variants detected in a cohort of 142 patients with premature ovarian insufficiency, diminished ovarian reserve, or unexplained infertility associated with normal ovarian function. We identified rare non-synonymous protein-altering variants in 2.8 % of women with ovarian deficiency and no such variants in our small cohort of women with infertility but normal ovarian function. We observed previously reported variants associated with premature ovarian insufficiency in patients with diminished ovarian reserve, highlighting a genetic relationship between these conditions. We confirmed functional impairment resulting from a p.Val15Met variant, detected for the first time in a patient with premature ovarian insufficiency. The remaining variants were associated with preserved transcriptional activity and localization of NR5A1, indicating that rare NR5A1 variants may be incorrectly curated if functional studies are not undertaken, and/or that NR5A1 variants may have only a subtle impact on protein function and/or confer risk of ovarian deficiency via oligogenic inheritance.

Author Correction: Human sex reversal is caused by duplication or deletion of core enhancers upstream of SOX9.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Human sex reversal is caused by duplication or deletion of core enhancers upstream of SOX9.

Disorders of sex development (DSDs) are conditions affecting development of the gonads or genitalia. Variants in two key genes, SRY and its target SOX9, are an established cause of 46,XY DSD, but the genetic basis of many DSDs remains unknown. SRY-mediated SOX9 upregulation in the early gonad is crucial for testis development, yet the regulatory elements underlying this have not been identified in humans. Here, we identified four DSD patients with overlapping duplications or deletions upstream of SOX9. Bioinformatic analysis identified three putative enhancers for SOX9 that responded to different combinations of testis-specific regulators. All three enhancers showed synergistic activity and together drive SOX9 in the testis. This is the first study to identify SOX9 enhancers that, when duplicated or deleted, result in 46,XX or 46,XY sex reversal, respectively. These enhancers provide a hitherto missing link by which SRY activates SOX9 in humans, and establish SOX9 enhancer mutations as a significant cause of DSD.

Mutant NR5A1/SF-1 in patients with disorders of sex development shows defective activation of the SOX9 TESCO enhancer.

Nuclear receptor subfamily 5 group A member 1/Steroidogenic factor 1 (NR5A1; SF-1; Ad4BP) mutations cause 46,XY disorders of sex development (DSD), with phenotypes ranging from developmentally mild (e.g., hypospadias) to severe (e.g., complete gonadal dysgenesis). The molecular mechanism underlying this spectrum is unclear. During sex determination, SF-1 regulates SOX9 (SRY [sex determining region Y]-box 9) expression. We hypothesized that SF-1 mutations in 46,XY DSD patients affect SOX9 expression via the Testis-specific Enhancer of Sox9 core element, TESCO. Our objective was to assess the ability of 20 SF-1 mutants found in 46,XY DSD patients to activate TESCO. Patient DNA was sequenced for SF-1 mutations and mutant SF-1 proteins were examined for transcriptional activity, protein expression, sub-cellular localization and in silico structural defects. Fifteen of the 20 mutants showed reduced SF-1 activation on TESCO, 11 with atypical sub-cellular localization. Fourteen SF-1 mutants were predicted in silico to alter DNA, ligand or cofactor interactions. Our study may implicate aberrant SF-1-mediated transcriptional regulation of SOX9 in 46,XY DSDs.

Comparative analysis of activins A and B in the adult mouse epididymis and vas deferens.

Activin A regulates testicular and epididymal development, but the role of activin B in the epididymis and vas deferens is unknown. Mouse models with reduced activin A (Inhba+/- and InhbaBK/+), or its complete absence (InhbaBK/BK), were investigated to identify specific roles of activins in the male reproductive tract. In 8-week-old Inhba+/- mice, serum activin A decreased by 70%, with a 50% reduction of gene expression and protein in the testis, epididymis and vas deferens. Activin B and the activin-binding protein, follistatin, were similar to wild-type. Testis weights were slightly reduced in Inhba+/- mice, but the epididymis and vas deferens were normal, while the mice were fertile. Activin A was decreased by 70% in the serum, testis, epididymis and vas deferens of InhbaBK/+ mice and was undetectable in InhbaBK/BK mice, but activin B and follistatin levels were similar to wild-type. In 6-week-old InhbaBK/BK mice, testis weights were 60% lower and epididymal weights were 50% lower than in either InhbaBK/+ or wild-type mice. The cauda epididymal epithelium showed infoldings and less intra-luminal sperm, similar to 3.5-week-old wild-type mice, but at 8 weeks, no structural differences in the testis or epididymis were noted between InhbaBK/BK and wild-type mice. Thus, Inhbb can compensate for Inhba in regulating epididymal morphology, although testis and epididymal maturation is delayed in mice lacking Inhba Crucially, reduction or absence of activin A, at least in the presence of normal activin B levels, does not lead to major defects in the adult epididymis or vas deferens.

Heat-induced masculinization in domesticated zebrafish is family-specific and yields a set of different gonadal transcriptomes.

Understanding environmental influences on sex ratios is important for the study of the evolution of sex-determining mechanisms and for evaluating the effects of global warming and chemical pollution. Fishes exhibit sexual plasticity, but the underlying mechanisms of environmental effects on their reproduction are unclear even in the well-established teleost research model, the zebrafish. Here we established the conditions to study the effects of elevated temperature on zebrafish sex. We showed that sex ratio response to elevated temperature is family-specific and typically leads to masculinization (female-to-male sex reversal), resulting in neomales. These results uncovered genotype-by-environment interactions that support a polygenic sex determination system in domesticated (laboratory) zebrafish. We found that some heat-treated fish had gene expression profiles similar to untreated controls of the same sex, indicating that they were resistant to thermal effects. Further, most neomales had gonadal transcriptomes similar to that of regular males. Strikingly, we discovered heat-treated females that displayed a normal ovarian phenotype but with a "male-like" gonadal transcriptome. Such major transcriptomic reprogramming with preserved organ structure has never been reported. Juveniles were also found to have a male-like transcriptome shortly after exposure to heat. These findings were validated by analyzing the expression of genes and signaling pathways associated with sex differentiation. Our results revealed a lasting thermal effect on zebrafish gonads, suggesting new avenues for detection of functional consequences of elevated temperature in natural fish populations in a global warming scenario.

Altered SOX9 genital tubercle enhancer region in hypospadias.

Human mutations in the SOX9 gene or its regulatory region can disrupt testicular development, leading to disorders of sex development (DSDs). Our previous work involving the genomic analysis of isolated DSD patients revealed a 78kb minimal sex determining region (RevSex) far upstream of SOX9 that was duplicated in 46,XX and deleted in 46,XY DSDs. It was postulated that RevSex contains a gonadal enhancer. However, the most highly conserved sub-region within RevSex, called SR4, was neither responsive to sex determining factors in vitro nor active in the gonads of transgenic mice, suggesting that SR4 may not be functioning as a testicular enhancer. Interestingly, SR4 transgenic mice showed reporter activity in the genital tubercle, the primordium of the penis and clitoris, a previously unreported domain of Sox9 expression. SOX9 protein was detected in the genital tubercle, notably in the urethral plate epithelium, preputial glands, ventral surface ectoderm and corpus cavernosa. SR4 may therefore function as a Sox9 genital tubercle enhancer, mutations of which could possibly lead to hypospadias, a birth defect seen in the DSD patients in the RevSex study. SR4 activity and the observed SOX9 expression pattern suggest that SR4 may function as a Sox9 genital tubercle enhancer. However, conditional ablation of Sox9 in the genital tubercle using Shh-Cre/+;Sox9flox/flox mice revealed no genital tubercle abnormalities, possibly due to compensation by similar Sox factors. To conclude, we have identified a novel regulatory enhancer driving Sox9 expression during external genitalia development.

Early depletion of primordial germ cells in zebrafish promotes testis formation.

As complete absence of germ cells leads to sterile males in zebrafish, we explored the relationship between primordial germ cell (PGC) number and sexual development. Our results revealed dimorphic proliferation of PGCs in the early zebrafish larvae, marking the beginning of sexual differentiation. We applied morpholino-based gene knockdown and cell transplantation strategies to demonstrate that a threshold number of PGCs is required for the stability of ovarian fate. Using histology and transcriptomic analyses, we determined that zebrafish gonads are in a meiotic ovarian stage at 14 days postfertilization and identified signaling pathways supporting meiotic oocyte differentiation and eventual female fate. The development of PGC-depleted gonads appears to be restrained and delayed, suggesting that PGC number may directly regulate the variability and length of gonadal transformation and testicular differentiation in zebrafish. We propose that gonadal transformation may function as a developmental buffering mechanism to ensure the reproductive outcome.